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1.  INTRODUCTION 


The  Naval  Biodynamics  Laboratory  (NBDL)  has  been  engaged  in  an  impact 
acceleration  injury  prevention  program  which  has  as  a  major  objective  the 
development  of  a  biodynamically  faithful  anthropomorphic  dummy.  It  is 
anticipated  that  this  development  will  provide  an  essential  component  to  any 
methodology  for  evaluating  proposed  restraint  and  protection  systems.  To  this 
end,  it  is  extremely  important  that  preliminary  research  determines  which 
variables  are  prime  correlates  with  injury  and  subsequently  establishes 
tolerance  limits,  in  terms  of  these  variables,  for  various  injury  levels. 

Since  1974,  NBDL  has  collected  an  extensive  data  base  from  -G^ 
accelerator  runs  on  Rhesus  and  human  subjects.  Because  of  safety 
considerations,  experiments  with  human  subjects  have  been  limited  to 
relatively  low-level  runs.  Experiments  with  Rhesus  subjects,  however,  have 
been  much  more  extensive,  and  some  have  resulted  in  fatal  injury. 

As  part  of  its  research  effort  in  this  area,  Desmatics,  Inc.  has  been 
actively  involved  in  the  development  of  statistically  based  models  to  predict 
injury  in  Rhesus  subjects  resulting  from  indirect  head/neck  impact 
acceleration  in  the  -X  direction.  Model  development  focused  initially  on 
fatality  as  the  primary  injury  criteria.  Subsequent  research,  however,  has 
emphasized  the  development  of  nonfatal  injury  prediction  models. 

Several  previous  Desmatics  technical  reports  ([2], [7], [8],  and  [9]) 
investigated  the  use  of  a  logistic  function  in  the  development  of  fatality 
prediction  models.  Empirical  data  was  used  to  fit  these  models  and  evaluate 
their  ability  to  accurately  predict  fatal  injuries.  The  variables  considered 
in  those  models  included  the  forces  and  torques  in  the  region  of  injury,  sled 
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profile  and  head  dynamic  response  variables,  and  variables  describing  the 
initial  position  of  the  head. 


The  data  base  analyzed  in  the  aforementioned  reports  consisted  of 
ninety-three  -G^  accelerator  runs  on  Rhesus  subjects,  conducted  between  1974 
and  1980  at  NBDL  [10].  During  1984  and  1985,  seventy-one  additional  -G^ 
accelerator  experiments  on  Rhesus  subjects  were  conducted  at  NBDL.  The 
objectives  of  the  present  technical  report  are: 

(1)  to  summarize,  update,  and  consolidate  the  Desmatics  modeling  efforts 
for  fatal  injury; 

(2)  to  investigate  extending  the  injury  prediction  models  to  include 
nonfatal  injury  criteria; 

and  (3)  to  evaluate  the  accuracy  and  validity  of  the  prediction  models 
obtained  in  the  modeling  efforts. 
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2.  METHODS 


The  logistic  model  has  the  following  functional  form:  J 

i 

_! 

P(x)  =  {1+expE-x'jJ]}  (2.1) 

where  x^  =  (1 ,Xj .x^ . . . ,x.  ) '  denotes  the  set  of  independent  (i.e.,  predictor) 
variables  considered, 

3_  =  (Sq,0j , . . . ,3^) '  denotes  the  set  of  parameter  values, 
and  P(x)  denotes  the  true  probability  of  injury  corresponding  to  jc. 

The  logistic  model  sometimes  is  motivated  by  threshold  theory.  For 
example,  suppose  that  each  Rhesus  subject  has  a  certain  unknown  threshold,  say 
T,  such  that  injury  will  occur  if  _>_  T  and  will  not  occur  if  x/ji  <  T.  If 
the  probability  distribution  of  Rhesus  thresholds  were  logistic  [which  is 
given  in  cumulative  distribution  form  in  equation  (2.1)],  model  (2.1)  would  be 
the  appropriate  model  according  to  threshold  theory.  It  may  be  noted  that  the 
logistic  distribution  closely  resembles  a  normal  distribution  but  has  proven 
to  be  more  analytically  tractable.  Furthermore,  the  assumption  that  the 
threshold  of  each  animal  can  reasonably  be  expressed  as  a  linear  combination 
of  predictor  variables  has  been  discussed  in  [6]. 

A  stepwise  selection  procedure  ([3]  and  [A]),  in  conjunction  with  model 
(2.1),  was  used  to  determine  the  inclusion  of  predictor  variables  in  model 
construction.  The  stepwise  selection  method  begins  with  no  variables  in  the 
model,  except  for  a  constant  term,  and  inserts  variables  in  turn  until  no 
further  improvement  can  be  made.  The  order  of  insertion  is  determined  by 
using  a  likelihood-ratio  test  [5]  as  a  measure  of  the  importance  of  variables 
not  yet  in  the  model  equation.  In  addition,  the  stepwise  procedure  can  also 


remove  a  variable  that  was  entered  at  an  earlier  stage  in  the  procedure. 
Accordingly,  judgements  on  the  insertion  or  removal  of  variables  are  made  at 
every  stage  in  the  stepwise  selection  procedure. 

The  models  presented  and  analyzed  in  this  report  were  based  on  two  sets 
of  injury  data.  The  first  of  these,  which  will  henceforth  be  referred  to  as 
Data  Set  A,  consisted  of  ninety-three  - G accelerator  runs  on  Rhesus  subjects, 
conducted  between  1974  and  1980  at  NBDL.  The  other  data  set,  which  will 
henceforth  be  referred  to  as  Data  Set  B,  consisted  of  seventy-one  -G^ 
accelerator  runs  on  Rhesus  subjects,  conducted  during  1984  and  1985  at  NBDL. 
Data  Sets  A  and  B  are  listed  in  Appendices  A  and  B,  respectively.  Table  1 
provides  a  brief  description  of  the  variables,  including  their  associated 
codes  and  units  of  measurement,  comprising  the  two  data  sets. 

In  Data  Set  A,  four  sled  profile  variables  (PSA,  ESV,  R00,  and  DOP),  two 
head  initial  condition  variables  (PHB  and  PHC),  and  nine  head  dynamic  response 
variables  (QHA,  QHB,  QHC,  RHA,  RHB,  RHC,  AAX,  AAY,  and  AAZ)  were  available  for 
analysis.  In  Data  Set  B,  these  fifteen  variables  and  an  additional  head 
initial  condition  variable  (PHA)  were  available  for  analysis.  Furthermore, 
because  of  symmetry  considerations,  absolute  values  of  several  of  the 
variables  were  used  in  model  development.  These  were  RHA,  RHC,  QHA,  QHC,  AAY, 
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PHA,  and  PHC. 

In  Data  Set  A,  two  sled  runs,  Run  Nos.  LX1363  and  LX1365,  were  omitted 
from  any  model  development.  Run  No.  LX1363,  which  was  a  fatal  experiment,  was 
omitted  due  to  suspicion  that  subject  animal  A03935  was  compromised  in  a 
previous  high-level  run  (Run  No.  LX1362)  made  earlier  that  day.  Run  No. 
LX1365,  which  also  was  a  fatal  experiment,  was  omitted  due  to  suspect  data 
quality.  Of  the  remaining  ninety-one  sled  experiments  in  Data  Set  A,  eleven 
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VARIABLE 


DESCRIPTION 


RUN  RUN  ID  NUMBER 

SUBJECT  ANIMAL  ID  NUMBER 

FATALITY  FATAL  (1)  OR  NONFATAL  (0) 

INJURY*  GRADED  INJURY  (0=No  injury,  l=Least  severe,  ....  4=Fatal) 

PSA  PEAK  SLED  ACCELERATION  (m/sec2) 

ESV  END  STROKE  VELOCITY  (m/sec) 

ROO  RATE  OF  SLED  ACCELERATION  ONSET  (m/sec3) 

DOP  DURATION  OF  PEAK  SLED  ACCELERATION  (msec) 

QHA*  PEAK  HEAD  ANGULAR  ACCELERATION,  X  COMPONENT  (rad/sec2) 

QHB  PEAK  HEAD  ANGULAR  ACCELERATION,  Y  COMPONENT  (rad/sec2) 

QHC*  PEAK  HEAD  ANGULAR  ACCELERATION,  Z  COMPONENT  (rad/sec2) 

RHA*  PEAK  HEAD  ANGULAR  VELOCITY,  X  COMPONENT  (rad/sec) 

RHB  PEAK  HEAD  ANGULAR  VELOCITY,  Y  COMPONENT  (rad/sec) 

RHC*  PEAK  HEAD  ANGULAR  VELOCITY,  Z  COMPONENT  (rad/sec) 

AAX  PEAK  HEAD  LINEAR  ACCELERATION,  X  COMPONENT  (m/sec2) 

AAY*  PEAK  HEAD  LINEAR  ACCELERATION,  Y  COMPONENT  (m/sec2) 

AAZ  PEAK  HEAD  LINEAR  ACCELERATION,  Z  COMPONENT  (m/sec2) 

PHA*+  INITIAL  ROLL  ANGLE  OF  THE  HEAD  (rad) 

PHB  INITIAL  PITCH  ANGLE  OF  THE  HEAD  (rad) 

PHC*  INITIAL  YAW  ANGLE  OF  THE  HEAD  (rad) 

♦absolute  value  of  variable  was  used  in  model  development 
variable  was  not  available  in  Data  Set  A 


Table  1.  Description  of  Variables. 


resulted  in  fatal  injury. 

In  Data  Set  B,  only  one  sled  run,  Run  No.  LX4824,  was  omitted  from  any 
model  development.  This  sled  experiment,  which  was  nonfatal,  was  deleted  due 
to  suspect  data  quality  resulting  from  sled  breakup  upon  impact.  Of  the 
remaining  seventy  sled  experiments  in  Data  Set  B,  only  two  resulted  in 
fatality. 


3.  FATALITY  MODELS 


This  section  summarizes  and  discusses  the  Desmatics  modeling  efforts 
regarding  the  development  of  fatality  prediction  models.  The  initial  phase  of 
these  efforts  focused  on  developing  and  comparing  fatality  prediction  models 
using  the  "historical"  data  base  (i.e.,  Data  Set  A).  The  resulting  models 
were  subsequently  used  to  obtain  predicted  probabilities  of  fatality  for  the 
more  recent  series  of  sled  experiments  (i.e.,  Data  Set  B). 

Because  Data  Set  B  had  only  two  fatal  sled  experiments,  there  was 
insufficient  information  to  attempt  to  develop  fatality  models  based  solely  on 
this  data  set.  Consequently,  the  final  phase  of  the  modeling  efforts  focused 
on  developing  fatality  models  using  all  available  data  (i.e.,  Data  Sets  A  and 
B  combined ) . 

3.1  Phase  I:  Models  Based  on  Data  Set  A 


Data  Set  A  contained  ninety-one  sled  experiments  and  fifteen  predictor 
variables  for  analysis.  Four  predictive  models  were  obtained  using  this  data 
set.  The  first  of  these  was  developed  from  the  set  of  four  sled  profile 
variables,  the  second  from  the  set  of  four  sled  profile  and  two  head  initial 
condition  variables,  the  third  from  the  set  of  nine  head  dynamic  response 
variables,  and  the  fourth  from  the  combined  set  of  fifteen  variables. 

It  should  be  pointed  out  that  different  sets  of  sled  experiments  were 
used  to  develop  each  of  these  four  models.  This  was  necessary  because  not 
every  predictor  variable  was  measured  in  each  sled  experiment.  For  example, 
there  were  a  number  of  sled  experiments  in  Data  Set  A  in  which  no  data  was 


available  on  any  of  the  head  dynamic  response  variables.  In  fact,  out  of  the 
ninety-one  sled  experiments ,  there  were  only  thirteen  in  which  the  complete 
set  of  fifteen  predictor  variables  was  measured.  In  general,  therefore,  the 
particular  set  of  sled  experiments  used  to  develop  a  prediction  model  depended 
on  the  variables  considered  for  inclusion  in  that  model. 

It  should  also  be  noted  that  some  of  the  Rhesus  subjects  were  run  more 
than  once.  Consequently,  some  dependencies  most  likely  exist  in  the  data  as  a 
result  of  having  multiple  sled  runs  with  the  same  animal.  However,  the 
assumption  is  made  here  that  any  resulting  bias  in  the  parameter  estimates  is 
small.  If  there  is  a  bias,  it  should  result  in  a  model  that  tends  to 
overpredict  probabilities. 

3.1.1  Sled  Profile  Variables 

Four  sled  profile  variables  were  considered  in  model  development:  PSA, 
ESV,  R00,  and  DOP.  It  was  found  that  the  complete  set  of  these  sled  variables 
was  available  on  eighty  sled  runs  in  Data  Set  A.  Of  these  eighty  runs,  ten 
were  fatal. 

The  "best"  one-variable  model,  as  determined  by  the  stepwise  selection 
procedure,  was  based  on  PSA,  peak  sled  acceleration.  There  was  no 
two-variable  model  involving  the  four  sled  profile  variables  which  offered  any 
significant  improvement  over  the  PSA-based  model. 

As  already  noted,  there  were  eighty  sled  experiments  in  Data  Set  A  in 
which  all  four  of  the  sled  profile  variables  were  measured.  PSA  alone, 
however,  was  available  for  every  sled  experiment  in  Data  Set  A.  Consequently, 
in  order  to  utilize  the  maximum  amount  of  information  available,  the  PSA-based 


model  was  subsequently  fit  to  the  complete  set  of  sled  experiments  in  Data  Set 
A.  That  model  was  given  by 

Pp  =  {l+exp[13.777-.012169(PSA)] }-1  (3.1) 

where  Pp  denotes  the  predicted  probability  of  fatal  injury. 

This,  in  fact,  illustrates  the  basic  approach  followed  in  obtaining  the 
models  presented  in  this  report.  First,  a  set  of  predictor  variables  is 
considered  for  inclusion  in  a  model.  Then,  the  stepwise  selection  procedure 
is  used  to  determine  which  variables  are  actually  to  be  included  in  the  final 
model.  The  selected  variables  are  then  fit  to  the  maximum  set  of  sled 
experiments  to  obtain  the  estimated  model  parameters. 

3.1.2  Sled  Profile  and  Head  Initial  Condition  Variables 

In  addition  to  the  four  sled  profile  variables,  two  head  initial 
condition  variables,  PHB  and  | PHC | ,  were  also  considered  in  model  development. 
There  were  seventeen  sled  runs,  which  included  ten  fatal  experiments, 
available  for  analysis  on  these  six  predictor  variables.  PHA,  another  head 
initial  condition  variable,  was  unavailable  in  Data  Set  A,  and  therefore  was 
not  included  in  this  analysis. 

As  determined  by  the  stepwise  selection  procedure,  the  best  one-variable 
model  was  based  on  PSA.  This  model  could  be  improved  on,  however,  by  the 
addition  of  ] PHC | ,  the  absolute  initial  yaw  angle  of  the  head.  Any  further 
additions  to  this  model  did  little  to  enhance  predictive  ability. 

The  two-variable  model  based  on  PSA  and  j PHC [  was  given  by 


Pp  =  {l+exp[17. 679-. 014587(PSA)-2. 5612 |PHC  j ] } 


(3.2) 


A 

where  Pp  denotes  the  predicted  probability  of  fatal  injury.  The  coefficients 
in  model  (3.2)  were  obtained  by  fitting  the  PSA  and  |PHC|  variables  to 
nineteen  sled  runs,  including  eleven  fatalities,  which  had  data  available  on 
those  two  variables. 

3.1.3  Head  Dynamic  Response  Variables 

Nine  head  dynamic  response  variables  were  considered  in  model 
development:  j QHA | ,  QHB,  |QHC| ,  | RHA | .  RHB,  |RHC|,  AAX,  [AAY | ,  and  AAZ.  A 

total  of  sixty-three  sled  runs,  which  included  seven  fatalities,  were 
available  for  analysis  on  these  nine  variables. 

The  best  one-variable  model,  as  determined  by  the  stepwise  selection 
procedure,  was  based  on  AAZ,  the  peak  Z-component  of  head  linear  acceleration. 
No  further  additions  resulted  in  any  significant  improvement  to  this  model. 

In  addition  to  the  aforementioned  sixty-three  runs,  one  additional  run 
(nonfatal)  had  data  available  for  the  AAZ  variable.  These  sixty-four  runs 
were  subsequently  used  to  obtain  the  following  AAZ-based  model: 

PF  =  {l+exp[7.03+.0053089(AAZ)]  }_1  (3.3) 

/\ 

where  Pp  denotes  the  predicted  probability  of  fatal  injury. 


3.1. A  Combined  Sled,  Head  Initial  Condition,  and  Dynamic  Response  Variables 


development.  There  were  thirteen  sled  runs,  which  included  seven  fatalities, 
available  for  analysis  on  the  combined  set  of  variables. 

The  best  one-variable  model,  as  determined  by  the  stepwise  selection 
procedure,  was  based  on  AAZ.  A  significantly  improved  model  was  subsequently 
achieved  by  the  addition  of  |PHCj  to  the  model.  Because  only  thirteen  sled 
experiments  were  used  in  the  analysis ,  enlarging  the  model  to  include  more 
than  two  variables  would  most  likely  result  in  an  overfit  model.  Such  models 
generally  fit  the  data  very  well  but  are  extremely  dangerous  for  prediction. 

Other  than  the  indicated  thirteen  sled  experiments,  no  additional  runs 
had  data  available  on  both  AAZ  and  PHC.  The  resulting  model  based  on  AAZ  and 
| PHC |  was  given  by 

PF  =  {l+exp[9.1391+.0056925(AAZ)-4.268|PHC|]}_1  (3.4) 

where  Pp  denotes  the  predicted  probability  of  fatal  injury. 

3.1.5  Comparison  of  Models 

The  four  fatality  prediction  models  presented  in  sections  3.1.1  through 
3.1.4  were  based  on  the  following  sets  of  predictor  variables: 


Model 

Variable ( s ) 

3.1 

PSA 

3.2 

PSA  and  | PHC | 

3.3 

AAZ 

3.4 

AAZ  and  1 PHC | 

where  PSA  denotes  the  peak  sled  acceleration,  PHC  denotes  the  initial  yaw 
angle  of  the  head,  and  AAZ  denotes  the  peak  Z-component  of  head  linear 


acceleration.  Unfortunately,  because  these  models  were  developed  from 
different  sets  of  sled  experiments,  they  cannot  be  compared  directly  with  each 
other.  For  example,  the  head  dynamic  response  model,  (3.3),  was  based  on  a 
set  of  sixty-four  sled  experiments.  However,  many  of  these  experiments  were 
relatively  low-level  (<40G)  runs,  which  are  inherently  much  easier  to  predict 
than  are  high-level  runs.  In  contrast,  the  thirteen  sled  experiments  used  to 
develop  the  combined-variables  model,  (3.4),  were  all  high-level  (>85G)  runs. 

It  was  noted,  however,  that  the  set  of  thirteen  sled  experiments  used  to 
develop  model  (3.4)  was  common  to  each  of  the  other  three  sets  of  sled 
experiments  used  in  model  development.  Thus,  some  assessment  of  the  relative 
performance  of  the  various  sets  of  predictor  variables  may  be  obtained  by 
comparing  the  predicted  probabilities  of  fatal  injury  for  these  thirteen  sled 
experiments.  Of  course,  in  order  to  ensure  a  common  basis  of  comparison,  the 
variables  in  models  (3.1),  (3.2),  and  (3.3)  were  refitted  specifically  to  this 
set  of  data. 

Performance  was  measured  using  the  Brier  statistic  [1]  which  is  defined 
as 

1  n  ^  2 
B  =  -  Z  (O.-P. ) 
n  .  ,  ii 
i=l 

t  h 

where  0^  is  the  observed  probability  of  fatality  (0  or  1)  for  the  i 
observation,  P^  is  the  corresponding  predicted  probability  of  fatality,  and  n 
is  the  number  of  observations.  The  Brier  statistic  expresses  the  average 
quadratic  loss;  thus,  a  smaller  value  of  B  indicates  better  performance.  In 
addition,  the  square  root  of  the  Brier  statistic  is  a  measure  of  the  average 
misprediction . 
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The  Brier  statistics  calculated  for  each  variable  set  were  as  follows: 


Set 

Variable (s) 

B 

ft 

#1 

PSA 

.212 

.460 

#2 

PSA  and  | PHC | 

.134 

.366 

#3 

AAZ 

.201 

.448 

#4 

AAZ  and  ] PHC | 

.148 

.385 

Not  surprisingly,  a  comparison  of  these  statistics  indicates  that  the  two 
two-variable  sets,  #2  and  #4,  were  better  than  either  of  the  two  one-variable 
sets,  #1  and  #3.  Furthermore,  variable  sets  #2  and  #4  were  fairly  comparable 
with  each  other,  as  were  sets  #1  and  #3. 

To  provide  a  "baseline"  for  comparison,  perhaps  the  simplest  prediction 
model  is  that  which  predicts  a  common  probability  of  fatality,  as  given  by  the 
observed  proportion  of  fatalities  in  the  data.  Of  the  thirteen  sled 
experiments,  53. 8Z  (i.e.,  7/13)  resulted  in  fatal  injury.  The  prediction 

A 

model  whereby  Pp=53.8Z  for  each  sled  experiment  has  an  associated  Brier 
statistic  of  B=.249.  As  indicated  above,  variable  set  #2  had  the  smallest 
Brier  statistic  given  by  B=.134.  This  represents  about  a  46Z  reduction  of  the 
baseline  value.  The  reduction  in  the  baseline  value  can,  in  some  sense,  be 
regarded  as  the  percentage  of  the  total  variation  in  the  data  that  can  be 
explained  by  the  predictor  variables.  Note  that  the  Brier  statistics  of  sets 
#1  and  #3,  which  are  the  two  one-variable  sets,  represent  only  a  15Z  and  a  19Z 
reduction,  respectively,  in  the  baseline  value.  The  use  of  variable  set  #4 
resulted  in  a  41Z  reduction  in  the  baseline  value. 

Regardless  of  which  variable  set  is  used,  a  large  percentage  of  the  total 
variation  in  the  data  remains  unexplained.  However,  it  should  be  noted  that 
each  of  the  thirteen  sled  experiments  under  consideration  was  close  to  the 


threshold  level  for  fatality.  These  runs  are  inherently  difficult  to  predict, 
especially  if  there  are  wide  differences  in  susceptibility  between  subjects. 
Since  there  is  no  data  available  for  estimating  differences  in  susceptibility, 
it  is  impossible  to  determine  at  this  time  whether  the  failure  to  more  fully 
explain  the  data  is  caused  by  inadequacy  of  the  models  or  by  intrinsic 
variability  resulting  from  differences  between  subjects. 

If  there  are  wide  differences  in  susceptibility  between  subjects,  those 
differences  must  be  taken  into  account  in  order  to  construct  a  good  prediction 
model.  Since  there  is  no  way  to  measure  susceptibility  directly,  it  must  be 
incorporated  into  the  model  indirectly  through  the  use  of  related  variables 
(e.g.,  anthropometric  variables).  Further  research  is  needed  in  order  to 
identify  such  variables. 

3.2  Predicted  Fatality  for  Data  Set  B 

The  four  fatality  prediction  models  (3.1),  (3.2),  (3.3),  and  (3.4)  were 
used  to  rank  the  sled  experiments  in  Data  Set  B  in  order  of  predicted 
severity.  These  rankings  are  provided  in  Tables  2a,  2b,  2c,  and  2d.  Included 
in  these  tables  are  the  lower  and  upper  95Z  confidence  limits  associated  with 
each  predicted  probability  of  fatality.  Also  included  is  the  level  of  graded 
injury  in  each  sled  experiment,  where  0=No  Injury,  l=Least  Severe,  .... 
4=Fatal.  Graded  injury  data  was  not  generally  available  for  the  sled 
experiments  in  Data  Set  A. 

There  were  only  two  fatalities  in  Data  Set  B.  These  were  Run  Nos.  LX4801 
and  LX4803.  As  can  be  seen  from  Tables  2a,  2b,  2c,  and  2d,  all  four  models 


failed  to  predict  that  LX4803  would  result  in  fatality.  On  the  other  hand, 


Lower  95X 

Predicted 

Upper  95X 

Graded 

Conf idence 

Probability 

Confidence 

Run 

Inlury 

PSA 

Limit 

of  Fatality 

Limit 

LX4787 

0 

90.1 

.00 

.00 

.01 

LX4792 

0 

95.5 

.00 

.00 

.01 

LX4819 

0 

100.5 

.00 

.00 

.01 

LX4809 

0 

100.6 

.00 

.00 

.01 

LX4813 

0 

100.7 

.00 

.00 

.01 

LX4800 

0 

100.9 

.00 

.00 

.01 

LX4821 

0 

100.9 

.00 

.00 

.01 

LX4815 

0 

101.0 

.00 

.00 

.01 

LX4807 

0 

101.1 

.00 

.00 

.01 

LX4823 

0 

101.1 

.00 

.00 

.01 

LX4811 

0 

101.2 

.00 

.00 

.01 

LX4802 

0 

101.6 

.00 

.00 

.01 

LX4817 

0 

101.8 

.00 

.00 

.01 

LX4804 

0 

102.0 

.00 

.00 

.01 

LX4798 

0 

104.3 

.00 

.00 

.01 

LX4782 

0 

124.2 

.00 

.00 

.01 

LX4781 

0 

124.5 

.00 

.00 

.01 

LX5122 

0 

197.1 

.00 

.00 

.02 

LX5159 

0 

199.4 

.00 

.00 

.02 

LX5127 

0 

200.1 

.00 

.00 

.02 

LX5137 

0 

200.5 

.00 

.00 

.02 

LX5134 

0 

201.6 

.00 

.00 

.02 

LX5140 

0 

202.7 

.00 

.00 

.02 

LX5149 

0 

202.9 

.00 

.00 

.02 

LX5143 

0 

203.4 

.00 

.00 

.02 

LX5146 

0 

203.6 

.00 

.00 

.02 

LX5154 

0 

204.6 

.00 

.00 

.02 

LX5123 

0 

403.7 

.00 

.00 

.05 

LX5128 

0 

406.4 

.00 

.00 

.05 

LX5132 

0 

408.2 

.00 

.00 

.05 

LX5135 

0 

409.2 

.00 

.00 

.05 

LX5151 

0 

410.8 

.00 

.00 

.05 

LX5150 

0 

411.0 

.00 

.00 

.05 

LX5157 

0 

411.1 

.00 

.00 

.05 

LX5155 

0 

412.0 

.00 

.00 

.05 

LX5152 

2 

413.5 

.00 

.00 

.05 

LX5138 

0 

413.9 

.00 

.00 

.05 

LX5156 

0 

415.2 

.00 

.00 

.05 

LX4808 

0 

544.4 

.00 

.00 

.11 

LX4820 

0 

545.0 

.00 

.00 

.11 

LX4810 

0 

545.2 

.00 

.00 

.11 

LX4788 

0 

552.2 

.00 

.00 

.11 

LX4783 

0 

556.0 

.00 

.00 

.11 

LX4784 

0 

625.1 

.00 

.00 

.17 

LX4814 

0 

692.5 

.00 

.00 

.25 

LX4789 

0 

699.0 

.00 

.01 

.26 

LX4818 

0 

700.0 

.00 

.01 

.26 

LX4816 

0 

700.3 

.00 

.01 

.26 

LX5125 

0 

720.0 

.00 

.01 

.30 

LX5160 

0 

723.9 

.00 

.01 

.30 

LX5129 

0 

724.5 

.00 

.01 

.30 

LX5131 

0 

724.5 

.00 

.01 

.30 

LX5147 

3 

727.7 

.00 

.01 

.31 

LX5141 

0 

730.6 

.00 

.01 

.31 

LX5144 

3 

731.6 

.00 

.01 

.32 

LX5164 

0 

732.5 

.00 

.01 

.32 

LX5161 

0 

733.0 

.00 

.01 

.32 

LX5162 

0 

734.4 

.00 

.01 

.32 

LX5165 

0 

744.3 

.00 

.01 

.34 

LX4785 

0 

777.1 

.00 

.01 

.41 

LX4790 

0 

834.4 

.00 

.03 

.55 

LX4806 

0 

838.5 

.00 

.03 

.56 

LX4822 

0 

842.6 

.00 

.03 

.57 

LX4812 

1 

845.4 

.00 

.03 

.58 

LX4803 

4 

884.1 

.00 

.05 

.67 

LX4786 

0 

940.0 

.00 

.09 

.80 

LX4795 

0 

985.1 

.00 

.14 

.88 

LX4791 

0 

988.4 

.00 

.15 

.88 

LX4801 

4 

994.7 

.00 

.16 

.89 

LX4799 

1 

1039.6 

.01 

.24 

.94 

Table  2a.  Fatality  Probabilities  in  Order  of  Predicted 
Severity  Based  on  Model  (3.1) 


Lower  95% 

Predicted 

Upper  95X 

Graded 

Confidence 

Probability 

Confidence 

Run 

Injury 

PSA 

IphcI 

Limit 

of  Fatality 

Limit 

LX4815 

0 

101.0 

• 

• 

• 

• 

LX4816 

0 

700.3 

• 

e 

e 

• 

LX4813 

0 

100.7 

.01 

.00 

.00 

.68 

LX4821 

0 

100.9 

.10 

.00 

.00 

.68 

LX4823 

0 

101.1 

.10 

.00 

.00 

.68 

LX4807 

0 

101.1 

.11 

.00 

.00 

.68 

LX4800 

0 

100.9 

.17 

.00 

.00 

.67 

LX4787 

0 

90.1 

.38 

.00 

.00 

.66 

LX4782 

0 

124.2 

.22 

.00 

.00 

.67 

LX4811 

0 

101.2 

.53 

.00 

.00 

.66 

LX4804 

0 

102.0 

.59 

.00 

.00 

.66 

LX5154 

0 

20 4.6 

.02 

.00 

.00 

.67 

LX3122 

0 

197.1 

.11 

.00 

.00 

.66 

LX4802 

0 

101.6 

.67 

.00 

.00 

.65 

LX5127 

0 

200.1 

.11 

.00 

.00 

.66 

LX5134 

0 

201.6 

.13 

.00 

.00 

.66 

LX4781 

0 

124.5 

.60 

.00 

.00 

.65 

LX5146 

0 

203.6 

.18 

.00 

.00 

.66 

LX4809 

0 

100.6 

.78 

.00 

.00 

.65 

LX3140 

0 

202.7 

.22 

.00 

.00 

.65 

LX5159 

0 

199.4 

.25 

.00 

.00 

.65 

LX5149 

0 

202.9 

.30 

.00 

.00 

.65 

LX5143 

0 

203.4 

.30 

.00 

.00 

.65 

LX5137 

0 

200.5 

.35 

.00 

.00 

.65 

LX4819 

0 

100.5 

.95 

.00 

.00 

.65 

LX4817 

0 

101.8 

1.06 

.00 

.00 

.66 

LX4798 

0 

104.3 

1.10 

.00 

.00 

.66 

LX4792 

0 

95.5 

1.60 

.00 

.00 

.69 

LX5135 

0 

409.2 

.02 

.00 

.00 

.63 

LX3123 

0 

403.7 

.30 

.00 

.00 

.62 

LX5128 

0 

406.4 

.32 

.00 

.00 

.61 

LX5157 

0 

411.1 

.30 

.00 

.00 

.62 

LX5155 

0 

412.0 

.32 

.00 

.00 

.61 

LX5132 

0 

408.2 

.40 

.00 

.00 

.61 

LX5156 

0 

415.2 

.38 

.00 

.00 

.61 

LX5152 

2 

413.5 

.45 

.00 

.00 

.61 

LX5138 

0 

413.9 

.48 

.00 

.00 

.61 

LX5151 

0 

410.8 

.50 

.00 

.00 

.61 

LX5150 

0 

411.0 

.52 

.00 

.00 

.61 

LX4783 

0 

556.0 

.00 

.00 

.00 

.62 

LX4820 

0 

545.0 

.21 

.00 

.00 

.60 

LX4788 

0 

552.2 

.30 

.00 

.00 

.59 

LX4808 

0 

544.4 

.45 

.00 

.00 

.59 

LX4810 

0 

545.2 

.75 

.00 

.00 

.58 

LX4784 

0 

625.1 

.30 

.00 

.00 

.58 

LX4814 

0 

692.5 

.18 

.00 

.00 

.58 

LX5144 

3 

731.6 

.00 

.00 

.00 

.60 

LX5164 

0 

732.5 

.00 

.00 

.00 

.60 

LX5125 

0 

720.0 

.08 

.00 

.00 

.59 

LX5160 

0 

723.9 

.07 

.00 

.00 

.59 

LX5165 

0 

744.3 

.00 

.00 

.00 

.59 

LX5162 

0 

734.4 

.06 

.00 

.00 

.59 

LX5161 

0 

733.0 

.07 

.00 

.00 

.59 

LX5129 

0 

724.5 

.14 

.00 

.00 

.58 

LX5141 

0 

730.6 

.15 

.00 

.00 

.58 

LX5131 

0 

724.5 

.20 

.00 

.00 

.58 

LX5147 

3 

727.7 

.30 

.00 

.00 

.57 

LX4785 

0 

777.1 

.08 

.00 

.00 

.58 

LX4818 

0 

700.0 

.65 

.00 

.00 

.56 

LX4789 

0 

699.0 

.80 

.00 

.00 

.57 

LX4806 

0 

838.5 

.15 

.00 

.01 

.57 

LX4822 

0 

842.6 

.15 

.00 

.01 

.57 

LX4803 

4 

884.1 

.10 

.00 

.01 

.58 

LX4790 

0 

834  .4 

.45 

.00 

.01 

.55 

LX4812 

1 

845.4 

.40 

.00 

.01 

.55 

LX4786 

0 

940.0 

.00 

.00 

.02 

.59 

LX4801 

4 

994.7 

.00 

.00 

.04 

.60 

LX4791 

0 

988.4 

.25 

.00 

.07 

.57 

LX4795 

0 

985.1 

.90 

.06 

.27 

.68 

LX4799 

1 

1039.6 

.80 

.13 

.39 

.73 

Table  2b.  Fatality  Probabilities  in  Order  of  Predicted  Severity 
Based  on  Model  (3.2) 


Lower  95X 

Predicted 

Upper  95X 

Graded 

Conf idence 

Probability 

Confidence 

Run 

Injury 

AAZ 

Limit 

of  Fatality 

Limit 

LX4815 

0 

• 

• 

• 

• 

LX4816 

0 

• 

• 

• 

• 

LX4786 

0 

• 

• 

• 

• 

LX4806 

0 

900. 

.00 

.00 

.01 

LX4795 

0 

700. 

.00 

.00 

.01 

LX4789 

0 

470. 

.00 

.00 

.02 

LX4788 

0 

400. 

.00 

.00 

.02 

LX4810 

0 

240. 

.00 

.00 

.03 

LX 4 8 08 

0 

200. 

.00 

.00 

.04 

LX4782 

0 

65. 

.00 

.00 

.05 

LX4807 

0 

30. 

.00 

.00 

.05 

LX4800 

0 

30. 

.00 

.00 

.05 

LX4792 

0 

28. 

.00 

.00 

.05 

LX4819 

0 

24. 

.00 

.00 

.05 

LX4811 

0 

21. 

.00 

.00 

.05 

LX4813 

0 

-18. 

.00 

.00 

.05 

LX4809 

0 

-27. 

.00 

.00 

.05 

LX4821 

0 

-32. 

.00 

.00 

.05 

LX4802 

0 

-35. 

.00 

.00 

.05 

LX5127 

0 

-42. 

.00 

.00 

.06 

LX4781 

0 

-43. 

.00 

.00 

.06 

LX4817 

0 

-43. 

.00 

.00 

.06 

LX4787 

0 

-52. 

.00 

.00 

.06 

LX4804 

0 

-57. 

.00 

.00 

.06 

LX4798 

0 

-60. 

.00 

.00 

.06 

LX4823 

0 

-70. 

.00 

.00 

.06 

LX5122 

0 

-70. 

.00 

.00 

.06 

LX5154 

0 

-82. 

.00 

.00 

.06 

LX5134 

0 

-82. 

.00 

.00 

.06 

LX5146 

0 

-85. 

.00 

.00 

.06 

LX5137 

0 

-92. 

.00 

.00 

.06 

LX5159 

0 

-100. 

.00 

.00 

.06 

LX5140 

0 

-125. 

.00 

.00 

.06 

LX5149 

0 

-140. 

.00 

.00 

.07 

LX5143 

0 

-180. 

.00 

.00 

.07 

LX5138 

0 

-190. 

.00 

.00 

.07 

LX5150 

0 

-205. 

.00 

.00 

.08 

LX5135 

0 

-215. 

.00 

.00 

.08 

LX5131 

0 

-220. 

.00 

.00 

.08 

LX5152 

2 

-240. 

.00 

.00 

.08 

LX5151 

0 

-240. 

.00 

.00 

.08 

LX5128 

0 

-280. 

.00 

.00 

.09 

LX5155 

0 

-280. 

.00 

• 

o 

o 

.09 

LX5157 

0 

-310. 

.00 

.00 

.09 

LX5156 

0 

-320. 

.00 

.00 

.09 

LX4783 

0 

-350. 

.00 

.01 

.10 

LX5123 

0 

-370. 

.00 

.01 

.10 

LX5132 

0 

-400. 

.00 

.01 

.11 

LX4822 

0 

-400. 

.00 

.01 

.11 

LX5129 

0 

-410. 

.00 

.01 

.11 

LX4820 

0 

-430. 

.00 

.01 

.12 

LX4814 

0 

-450. 

.00 

.01 

.12 

LX4818 

0 

-450. 

.00 

.01 

.12 

LX5164 

0 

-500. 

.00 

.01 

.13 

LX5160 

0 

-500. 

.00 

.01 

.13 

LX5141 

0 

-500. 

.00 

.01 

.13 

LX4790 

0 

-500. 

.00 

.01 

.13 

LX5162 

0 

-510. 

.00 

.01 

.13 

LX4799 

1 

-510. 

.00 

.01 

.13 

LX5161 

0 

-520. 

.00 

.01 

.14 

LX5125 

0 

-550. 

.00 

.02 

.15 

LX4803 

4 

-550. 

.00 

.02 

.15 

LX4791 

0 

-600. 

.00 

.02 

.16 

LX4812 

1 

-610. 

.00 

.02 

.16 

LX4784 

0 

-680. 

.00 

.03 

.19 

LX5144 

3 

-710. 

.01 

.04 

.20 

LX5147 

3 

-850. 

.02 

.07 

.27 

LX4785 

0 

-1100. 

.09 

.23 

.50 

LX5165 

0 

-1200. 

.14 

.34 

.63 

LX4801 

4 

-2250. 

.70 

.99 

1.00 

Table  2c.  Fatality  Probabilities  in  Order  of  Predicted 
Severity  Based  on  Modpl  (3.3) 


Lower  95%  Predicted  Upper  95% 
Graded  Confidence  Probability  Confidence 


Run 

Injury 

AAZ 

iPHCl 

Limit 

of  Fatality 

LX4815 

0 

• 

e 

e 

• 

• 

LX4816 

0 

e 

e 

e 

• 

• 

LX4786 

0 

• 

• 

e 

• 

• 

LX4806 

0 

900. 

.15 

.00 

.00 

.94 

LX4788 

0 

400. 

.30 

.00 

.00 

.89 

LX4795 

0 

700. 

.90 

.00 

.00 

.94 

LX4813 

0 

-18. 

.01 

.00 

.00 

.84 

LX4807 

0 

30. 

.11 

.00 

.00 

.84 

LX5154 

0 

-82. 

.02 

.00 

.00 

.83 

LX4800 

0 

30. 

.17 

.00 

.00 

.84 

LX4782 

0 

65. 

.22 

.00 

.00 

.85 

LX4821 

0 

-32. 

.10 

.00 

.00 

.83 

LX5127 

0 

-42. 

.11 

.00 

.00 

.83 

LX4789 

0 

470. 

.80 

.00 

.00 

.91 

LX4808 

0 

200. 

.45 

.00 

.00 

.87 

LX4823 

0 

-70. 

.10 

.00 

.00 

.83 

LX5122 

0 

-70. 

.11 

.00 

.00 

.83 

LX5134 

0 

-82. 

.13 

.00 

.00 

.82 

LX5146 

0 

-85. 

.18 

.00 

.00 

.82 

LX5135 

0 

-215. 

.02 

.00 

.00 

.80 

LX5159 

0 

-100. 

.25 

.00 

.00 

.82 

LX5140 

0 

-125. 

.22 

.00 

.00 

.81 

LX4810 

0 

240. 

.75 

.00 

.00 

.89 

LX4787 

0 

-52. 

.38 

.00 

.00 

.83 

LX4783 

0 

-350. 

.00 

.00 

.00 

.78 

LX5137 

0 

-92. 

.35 

.00 

.00 

.82 

LX5149 

0 

-140. 

.30 

.00 

.00 

.81 

LX5131 

0 

-220. 

.20 

.00 

.00 

.80 

LX4811 

0 

21. 

.53 

.00 

.00 

.84 

LX5143 

0 

-180. 

.30 

.00 

.00 

.80 

LX4781 

0 

-43. 

.60 

.00 

.00 

.84 

LX4804 

0 

-57. 

.59 

.00 

.00 

.83 

LX5164 

0 

-500. 

.00 

.00 

.00 

.75 

LX4822 

0 

-400. 

.15 

.00 

.00 

.76 

LX5129 

0 

-410. 

.14 

.00 

.00 

.76 

LX5128 

0 

-280. 

.32 

.00 

.00 

.78 

LX5155 

0 

-280. 

.32 

.00 

.00 

.78 

LX5157 

0 

-310. 

.30 

.00 

.00 

.78 

LX4802 

0 

-35. 

.67 

.00 

.00 

.84 

LX5138 

0 

-190. 

.48 

.00 

.00 

.80 

LX5160 

0 

-500. 

.07 

.00 

.00 

.75 

LX5162 

0 

-510. 

.06 

.00 

.00 

.75 

LX5161 

0 

-520. 

.07 

.00 

.00 

.74 

LX5152 

2 

-240. 

.45 

.00 

.00 

.79 

LX4814 

0 

-450. 

.18 

.00 

.00 

.75 

LX4820 

0 

-430. 

.21 

.00 

.00 

.75 

LX5150 

0 

-205. 

.52 

.00 

.00 

.80 

LX5123 

0 

-370. 

.30 

.00 

.00 

.76 

LX5156 

0 

-320. 

.38 

.00 

.00 

.78 

LX5125 

0 

-550. 

.08 

.00 

.00 

.74 

LX4809 

0 

-27. 

.78 

.00 

.00 

.85 

LX5141 

0 

-500. 

.15 

.00 

.00 

.74 

LX5151 

0 

-240. 

.50 

.00 

.00 

.80 

LX4803 

4 

-550. 

.10 

.00 

.00 

.74 

LX4819 

0 

24. 

.95 

.00 

.01 

.88 

LX5132 

0 

-400. 

.40 

.00 

.01 

.76 

LX5144 

3 

-710. 

.00 

.00 

.01 

.72 

LX4791 

0 

-600. 

.25 

.00 

.01 

.72 

LX4790 

0 

-500. 

.45 

.00 

.01 

.74 

LX4817 

0 

-43. 

1.06 

.00 

.01 

.89 

LX4798 

0 

-60. 

1.10 

.00 

.02 

.90 

LX4784 

0 

-680. 

.30 

.00 

.02 

.70 

LX4812 

1 

-610. 

.40 

.00 

.02 

.72 

LX4818 

0 

-450. 

.65 

.00 

.02 

.77 

LX5147 

3 

-850. 

.30 

.00 

.05 

.67 

LX4799 

1 

-510. 

.80 

.00 

.06 

.79 

LX4785 

0 

-1100. 

.08 

.00 

.07 

.67 

LX4792 

0 

28. 

1.60 

.00 

.08 

.97 

LX5165 

0 

-1200. 

.00 

.00 

.09 

.69 

LX4801 

4 

-2250. 

.00 

.11 

.98 

1.00 

Table  2d.  Fatality  Probabilities  in  Order  of  Predicted  Severity 
Based  on  Model  (3.41 


LX A 801  was  accurately  predicted  by  models  (3.3)  and  (3.4),  but  was  poorly 
predicted  by  models  (3.1)  and  (3.2).  It  may  also  be  noted  that  all  models  had 
fairly  wide  confidence  intervals,  as  expressed  by  the  difference  between  the 
upper  and  lower  confidence  limits.  Model  (3.3),  however,  provided  the 
narrowest  intervals. 

The  confidence  intervals  in  Tables  2a-2d  are  wide  for  two  basic  reasons. 
First,  dichotomous  variables  are  inherently  difficult  to  predict  since  so 
little  information  is  obtained  from  each  experiment.  In  general,  much  larger 
sample  sizes  are  required  for  these  variables  than  would  be  required  to  obtain 
equivalent  precision  for  continuous  variables.  The  second  problem  arises  from 
the  fact  that  much  of  the  data  used  to  build  these  models  was  obtained  from 
relatively  low-level  runs  and  thus  contributed  little  to  the  model-building 
process.  The  runs  which  are  most  useful  for  estimating  the  parameters  of  a 
logistic  model  are  those  which  are  fairly  close  to  the  fatality  threshold. 

The  Brier  statistics  for  the  four  models  were  as  follows: 


Model 

B 

(3.1) 

.025 

.157 

(3.2) 

.031 

.177 

(3.3) 

.017 

.131 

(3.4) 

.015 

.123 

The  baseline  Brier  value  was  B=.028.  As  indicated  from  the  statistics 
above,  models  (3.3)  and  (3.4)  were  clearly  better  than  either  of  (3.1)  or 
(3.2).  These  latter  two  models  resulted  in  little  or  no  reduction  in  the 
baseline  Brier  value.  Model  (3.4)  had  a  Brier  value  which  was  slightly  less 
than  that  for  model  (3.3).  It  should  be  expected,  however,  that  model  (3.4) 
would  outperform  model  (3.3)  since  (3.4)  contains  an  additional  predictor 


variable.  Furthermore,  the  Brier  statistic  of  B=.015  associated  with  model 
(3.4)  represented  a  46. 42  reduction  of  the  baseline  value.  This  reduction  is 
consistent  with  that  found  in  section  3.1.5. 

It  should  be  noted  that  the  superior  performance  of  the  two  models  which 
included  as  a  predictor  the  peak  Z-component  of  head  linear  acceleration  (3.3 
and  3.4)  can  to  a  large  extent  be  attributed  to  their  predictions  for  one  sled 
run.  That  run  (LX4801)  had  a  moderately  high  peak  sled  acceleration  but  an 
extremely  high  head  acceleration.  Extreme  caution  is  needed  when  attempting 
to  draw  conclusions  from  a  single  experimental  run. 

With  regard  to  nonfatal  injury,  it  is  not  unreasonable  to  expect  that  a 
fatality  model  should  predict  nonfatal  injury  to  some  extent,  assuming  that 
the  failure  mechanism  is  the  same  and  that  injury  is  more  or  less  a  matter  of 
degree.  It  is  difficult  to  state,  however,  what  magnitude  of  probability  that 
a  fatality  model  should  predict  for  a  specific  level  of  nonfatal  injury. 

As  can  be  seen  from  Tables  2a,  2b,  2c,  and  2d,  sled  experiments  with 
nonzero  graded  injury  levels  were  generally  ranked  higher  in  terms  of  their 
predicted  probability  of  fatality  than  were  sled  experiments  involving  no 
injury.  The  actual  ordering  of  the  injuries  was  not  entirely  consistent, 
however,  with  their  severity.  For  example,  fatalities  were  not  always  ranked 
higher  than  level-3  injuries,  which  were  not  always  ranked  higher  than  level-2 
injuries,  and  so  forth.  It  was  found,  though,  that  model  (3.3)  provided  the 
most  consistent  ranking  of  sled  experiments  relative  to  their  associated 
injury  severity. 


Desmatics  combined  Data  Sets  A  and  B  and  consolidated  its  modeling 
efforts.  As  was  done  previously  in  Section  3.1,  models  were  developed  from 
the  sled  profile  variables,  the  sled  profile  and  head  initial  condition 
variables,  the  dynamic  response  variables,  and  all  variables  combined. 

Although  PHA,  the  initial  roll  angle  of  the  head,  was  measured  in  Data  Set  B, 
it  was  omitted  here  because  it  was  not  available  for  sled  experiments  in  Data 
Set  A;  otherwise,  the  stepwise  selection  procedure  would  have  been  limited  to 
Data  Set  B  anytime  PHA  was  considered  for  inclusion  in  the  model. 

As  before,  the  stepwise  selection  method  was  used  to  determine  the 
inclusion  of  predictor  variables  in  the  various  models.  Except  for  the  model 
developed  from  the  set  of  sled  profile  and  head  initial  condition  variables, 
the  predictor  sets  which  were  selected  were  the  same  as  those  determined 
previously  in  Section  3.1.  In  the  case  of  the  sled  profile  and  head  initial 
condition  variables,  the  model  based  on  PSA  and  |PHC|  was  replaced  by  a  model 
based  on  PSA  and  PHB,  where  PHB  is  the  initial  pitch  angle  of  the  head. 

Maximum  closure  on  these  models  was  obtained  by  refitting  the  respective 
variable  sets  to  the  complete  set  of  sled  experiments  in  Data  Sets  A  and  B 
which  had  data  available  on  those  variables.  The  fitted  models  are  summarized 
in  Table  3.  For  the  sake  of  completeness  and  comparison,  the  model  based  on 
PSA  and  |PHC|  is  also  included  in  this  table. 

3.3.1  Comparison  of  Models 


Models  (3.5)  through  (3.9)  listed  in  Table  3  are  not  directly  comparable 


because  they  were  not  fit  to  the  same  set  of  sled  experiments.  However,  the 
set  of  eighty  sled  experiments  used  to  fit  model  (3.9)  was  common  to  the 
construction  of  each  of  the  other  four  models.  Thus,  in  order  to  provide  a 
common  basis  for  comparison,  models  (3.5)  through  (3.8)  were  refitted 


specifically  to  these  eighty  runs.  The  corresponding  Brier  statistics  were  as 
follows : 


Variable (s ) 

B 

Reduction  of 

PSA 

.054 

.232 

46% 

PSA  and  | PHC| 

.052 

.228 

48% 

PSA  and  PHB 

.047 

.217 

53% 

AAZ 

.049 

.221 

51% 

AAZ  and  1  PHCl 

.037 

.192 

63% 

The  baseline  value  of  the  Brier  statistic  was  B=.100. 

As  seen  from  the  statistics  above,  AAZ  was  slightly  better  than  PSA  as 
the  best  single  predictor  of  fatality.  In  addition,  the  AAZ-based  model  was 
found  to  have  narrower  confidence  limits  than  the  PSA-based  model.  The  best 
two-variable  model,  and  best  model  overall,  was  that  based  on  AAZ  and  |PHC|, 


This  model  resulted  in  a  63%  reduction  of  the  baseline  Brier  value. 
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3.3.2  Tolerance  Thresholds 


An  injury  tolerance  threshold  of  t%,  relative  to  a  given  prediction 
model,  is  defined  as  the  set  of  values  of  the  predictor  variables  where  injury 
is  predicted  for  t%  of  the  subjects,  i.e.,  where  the  predicted  probability  of 
injury  is  t%.  An  tolerance  threshold  of  t%  will  be  denoted  by  LDt.  In  the 
rase  of  a  fatality  prediction  model,  the  LD^  represents,  of  course,  a  fatality 
tolerance  threshold. 

Table  4  shows  the  LD^s  corresponding  to  models  (3.5)  through  (3.9)  for 
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Model 

(3.5) 

Model 

(3.6) 

Model 

(3.7) 

Model 

(3.8) 

Model 

(3.9) 

“5 

x1=854.7 

x2=8.421 

x3=9.083 

x . =-734 .0 

4 

x5=4.022 

LD10 

x1=922.2 

x2=9.168 

x3=9.830 

x, =-892.5 

4 

x5=4.769 

“>25 

Xj=1021.4 

x2=10.266 

x3=10.928 

x. =-1125. 5 

4 

x5=5.867 

LD50 

x1=1120.7 

x2=11.365 

x3=12.027 

x . =-1358 .7 

4 

x5=6.966 

LD75 

x1=1220.0 

x2=12.464 

x3=13.126 

x. =-1592.0 

4 

x3=8.065 

LD90 

x1=1319.2 

x2=13.562 

x3=14.224 

x. =-1825.0 

4 

x5=9.163 

LD95 

x1=1386.7 

x2=14.309 

x3=14.971 

x. =-1983. 5 

4 

x5=9.910 

KEY:  x1=PSA 

x2= .0098119 ( PSA) +.88201 |PHC | 

x3= .011217 ( PSA ) -2 . 1569 ( PHB ) 

x. =AAZ 
4 

x5=-.00497(AAZ)+2.7425 | PHC | 

Table  4.  Fatality  Tolerance  Thresholds 
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t=5, 10, 25, 50, 75, 90,  and  95.  It  should  be  noted  that  for  a  single-variable 

model,  the  LDt  can  be  calculated  uniquely.  For  example,  the  LD^  (l.e.,  50t 

2 

chance  of  fatality)  for  model  (3.5)  is  given  by  PSA=1120.7  m/sec  =  114.2  G. 
For  a  two-variable  model,  the  LDfc  is  represented  by  the  equation  of  a  straight 
line  in  the  respective  variables.  For  example,  the  LD^  for  model  (3.9)  is 
given  by  any  values  of  AAZ  and  |PHC|  which  satisfy  the  equation 
6.9663  =  -.00497(AAZ)  +  2.7425 | PHC| . 


4.  NON FATAL  INJURY  MODELS 


Because  graded  injury  information  was  not  generally  available  for  those 
sled  experiments  in  Data  Set  A,  the  development  of  nonfatal  injury  prediction 
models  was  necessarily  based  on  an  analysis  of  Data  Set  B.  There  was  the 
following  injury  information  in  that  data  set: 

Graded 

Run  No.  Injury 


LX4799  1 
LX4801  4 
LX4803  4 
LX4812  1 
LX5144  3 
LX5147  3 
LX5152  2 
Other  0 


where  0=No  Injury,  l=Least  Severe,  ...»  4=Fatal. 

The  use  of  the  logistic  model  requires  that  the  response  variable  has 
either  a  value  of  0  (i.e.,  no  injury  occurs)  or  1  (i.e.,  injury  occurs). 
Accordingly,  the  following  four  injury  response  variables  could  be  defined: 
To  if  the  level  of  graded  injury  <  j 


I.  =j 

J  (l  if  the  level  of  graded  injury  >_  j 
for  j=l,2,3,  and  4.  Thus,  for  example,  1^  =  1  if  any  (nonzero)  level  of  injury 
occurs  and  1^=0  if  no  injury  occurs.  Similarly,  1^=0  if  no  fatality  occurs 
and  1^=1  if  a  fatality  occurs.  Correspondingly,  let  P_.  denote  the  probability 
that  I  =1.  Equivalently,  represents  the  probability  that  a  graded  injury 
of  level  j  or  greater  will  occur. 

There  are  two  primary  issues  involving  the  development  of  nonfatal  injury 


prediction  models  which  need  to  be  addressed.  First,  the  prime  correlates  of 
nonfatal  injury  need  to  be  determined  and  compared  with  those  determined 


earlier  in  fatality  model  development.  Secondly,  the  ability  of  these  models 
to  "separate"  or  "distinguish"  between  the  various  injury  levels  must  be 
ascertained.  For  instance,  it  is  reasonable  to  require  that  the  derived 

A  A  A  A 

models  have  the  property  that  Nonfatal  injury  prediction  models 

would  be  of  little  use  if  separation  of  injury  levels  cannot  be  achieved. 

4.1  Model  Development 


The  development  and  construction  of  nonfatal  injury  prediction  models 
proceeded  according  to  the  approach  followed  in  Section  3  when  fatality  was 
the  dependent  variable  of  interest.  Using  1^ ,  I 2,  and  as  dependent 
variables,  prediction  models  were  developed  from  the  set  of  sled  profile 
variables,  sled  profile  and  head  initial  condition  variables,  head  dynamic 
response  variables,  and  all  variables  combined.  The  important  predictor 
variables,  as  determined  by  the  stepwise  selection  procedure,  were  essentially 
the  same  as  those  selected  previously  in  the  construction  of  fatality 
prediction  models.  This  result  was  encouraging  since  it  suggested  that  the 
same,  or  similar,  mechanisms  that  underlie  fatality  may  also  underlie  nonfatal 
injury. 

AAZ,  the  peak  Z-component  of  head  linear  acceleration,  was  determined 
earlier  to  be  the  best  single  predictor  of  fatality.  This  variable  was 
similarly  found  to  be  the  best  single  predictor  of  the  response  variables  1^, 
I^,  and  I^.  The  following  AAZ-based  nonfatal  injury  models  were  obtained: 
P^fl+expCS^TO+.OOSeOASCAAZM  f1,  (4.1) 

P2=(l+exp[4.0526+.0034072(AAZ)] f1 ,  (4.2) 

and  P3={l+exp[4.9206+.0042952(AAZ)]}~1,  (4.3) 
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AAZ  was  measured  on  sixty-seven  of  the  seventy  sled  experiments  in  Data  Set  B. 

Models  (4.1),  (4.2),  and  (4.3)  were  fit  to  these  sixty-seven  runs,  involving 

seven,  five,  and  four  injuries,  respectively.  For  predicting  fatality  (i.e., 

1^),  model  (3.8),  which  was  developed  in  Section  3,  may  be  used: 

P4=PF={l+exp[6.4031+.0047125(AAZ)3}_1.  (3.8) 

As  an  indication  of  whether  these  models  provided  separation  of  the 

various  injury  levels,  their  LD^s  may  be  compared.  The  LD,.qS  of  models 

(4.1),  (4.2),  (4.3),  and  (3.8)  are  -1011.7,  -1189.4,  -1145.6,  and  -1358.7 
2 

m/sec  ,  respectively.  Ideally,  the  magnitudes  of  these  tolerance  thresholds 
should  increase  as  the  injury  level  increases.  As  indicated,  however,  the 
magnitudes  are  not  strictly  increasing.  Specifically,  the  j  |  of  P^  is 
less  than  that  of  P2« 

Evidently,  a  graded  injury  of  level-3  severity  cannot  be  separated,  in  a 
predictive  sense,  from  an  injury  of  level-2  severity.  (This  is  not  surprising 
since  there  was  only  one  level-2  injury  in  the  data  base.)  Therefore,  these 
two  levels  of  injury,  at  least  for  this  data  set,  should  probably  be 
considered  to  be  of  equal  severity.  This  would  imply,  then,  that  until 
further  injury  information  becomes  available  and  can  be  analyzed,  model  (4.2) 
should  be  used  to  predict  both  level-2  and  level-3  injuries. 

Models  (4.1),  (4.2),  and  (3.8)  are  illustrated  graphically  in  Figure  1. 

As  can  be  seen  from  this  figure,  there  is  still  not  complete  separation  of  P^ , 

/N  A 

P2>  and  Pp,  primarily  at  high  negative  AAZ  levels.  For  the  most  part, 
however,  there  is  adequate  separation  for  prediction.  Note  that  the  crossover 
does  not  occur  until  the  predicted  probability  exceeds  90%. 

Continuing  with  its  analysis,  Desmatics  attempted  to  extend  the  AAZ-based 
injury  models  to  two-variable  models.  Using  1^ ,  I2,  and  as  the  injury 
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Figure  1.  Graphical  Illustration  of  P 


response  variables  it  was  found  that  no  two-variable  model  provided  any 
significant  improvement  over  the  predictive  abilities  of  the  AAZ-based  models. 
This  result  was  not  too  surprising  owing  to  the  rather  limited  injury 
information  available  in  Data  Set  B. 

4,2  Injury  Predictions  for  Data  Set  A 

As  previously  mentioned,  graded  injury  data  was  unavailable  for  Data  Set 
A.  However,  for  completeness,  models  (4.1),  (4.2),  and  (3.8)  were  used  to 
obtain  nonfatal  injury  predictions  for  the  ninety-one  sled  experiments  in  Data 
Set  A.  Table  5  shows  the  P  ,  P  ,  and  P_  predictions  which  were  obtained. 


Run  Fatality 
LXA5A  0 

LX650  0 

LX651  0 

LX653  0 

LX65A  0 

LX656  0 

LX 657  0 

LX659  0 

LX660  0 

LX661  1 

LX3006  0 

LX3007  0 

LX3008  0 

LX3009  0 

LX3010  0 

LX 30 11  0 

LX3012  0 

LX3013  0 

LX301A  0 

LX3015  0 

LX3016  1 

LX3017  0 

LX3018  0 

LX3019  0 

LX3020  0 

LX3031  1 

LX3033  1 


LX3027  0  87 
LX3185  0  62 
LX3025  0  A2 
LX3708  0  33 
LX370A  0  27 
LX3696  0  12 
LX3703  0  -13 
LX3713  0  -20 
LX3692  0  -26 
LX1081  0  -27 
LX3700  0  -30 
LX302A  0  -A 5 
LX3695  0  -60 
LX3187  0  -63 
LX3691  0  -67 
LX3189  0  -88 
LX3183  0  -90 
LX3699  0  -102 
LX3707  0  -105 
LX3191  0  -108 
LX3710  0  -117 
LX1899  0  -165 


AAZ  Pi  P2 


• 

.02 

.01 

.02 

.01 

.02 

.01 

.02 

.02 

.02 

.02 

.02 

.02 

.03 

.02 

.03 

.02 

.03 

.02 

.03 

.02 

.03 

.02 

.03 

.02 

.03 

.02 

.03 

.02 

.03 

.02 

.03 

.02 

.03 

.02 

.0A 

.02 

.OA 

.02 

.OA 

.02 

.OA 

.03 

.05 

.03 

Table  5.  Injury  Predictions  for  Data  Set  A  Sled  Experiments 


5.  SUMMARY  AND  DISCUSSION 


This  report  documents  the  construction  and  evaluation  of  a  set  of 
indirect  impact  injury  prediction  models  for  Rhesus  monkeys  under  -G^ 
acceleration.  Two  data  bases  were  used  in  model  construction:  (1)  a  set  of 
ninety-three  accelerator  runs  conducted  between  1974  and  1980  at  NBDL  (Data 
Set  A)  and  (2)  a  set  of  seventy-one  additional  runs  conducted  at  NBDL  during 
1984  and  1985  (Data  Set  B). 

Portions  of  Data  Set  A  have  been  analyzed  and  discussed  in  previous 
Desmatics  reports.  The  first  step  in  the  current  effort  was  to  summarize, 
consolidate,  and  update  those  earlier  modeling  efforts.  Four  fatality 
prediction  models  were  fit  to  the  data,  using  a  stepwise  variable  selection 
procedure.  The  first  model  was  developed  from  sled  profile  variables  only, 
the  second  from  sled  profile  and  head  initial  condition  variables,  the  third 
from  head  dynamic  response  variables,  and  the  fourth  from  the  combined  set  of 
all  of  the  variables. 

Peak  sled  acceleration  (PSA)  was  found  to  be  the  best  predictor  of 
fatality  from  among  the  sled  profile  variables,  and  the  inclusion  of 
additional  sled  variables  did  not  substantially  improve  the  prediction  model. 
Of  the  head  dynamic  response  variables,  the  peak  Z-component  of  head  linear 
acceleration  (AAZ)  was  the  best  predictor,  and  again  the  model  could  not  be 
substantially  improved  by  the  inclusion  of  additional  variables.  However,  it 
was  found  that  either  of  these  models  could  be  improved  significantly  by 
including  a  head  initial  condition  variable,  the  absolute  initial  yaw  angle  of 
the  head  ( | PHC | ) . 

Each  of  the  fatality  prediction  models  developed  from  Data  Set  A  was 


refit  to  the  set  of  thirteen  accelerator  runs  for  which  all  of  the  relevant 
predictor  variables  were  available.  Performance  of  the  models  was  evaluated 
using  Brier  statistics,  which  measure  the  average  squared  difference  between 
the  observations  (0  or  1)  and  the  predicted  probabilities.  The  two 
one-variable  models,  based  on  PSA  and  AAZ,  respectively,  were  found  to  be 
fairly  comparable,  as  were  the  two  two-variable  models,  which  included  |PHC|. 

One  measure  of  model  performance  is  the  percentage  reduction  in  the  Brier 
statistic  when  the  model  is  fit,  as  compared  to  a  baseline  value  obtained  by 
predicting  a  constant  average  probability  of  fatality  for  each  accelerator 
run.  This  reduction  may  be  thought  of  as  a  measure  of  the  amount  of 
variability  in  the  data  which  can  be  accounted  for  by  the  variables  in  the 
model.  By  this  criteria,  the  one-variable  models  were  able  to  explain  less 
than  20Z  of  the  variability  in  the  data,  and  the  two-variable  models  only 
about  50Z . 

It  may  not  be  possible  to  develop  better  prediction  models  using  only  the 
predictor  variables  now  available,  since  prediction  accuracy  is  limited  by  the 
inherent  variability  in  the  subject  population.  If  there  are  wide  differences 
in  susceptibility  between  subjects,  no  model  can  perform  well  unless  it 
incorporates  some  measure  of  those  differences.  In  any  case,  some  measure  of 
between-subject  variability  is  needed  before  the  current  models  can  be  truly 
evaluated.  One  way  to  obtain  such  an  estimate  would  be  to  test  several 
subjects  under  identical  conditions  and  compare  the  results.  It  might  be 
possible  to  relate  those  results  to  some  measurable  anthropometric  variables 
which  could  then  be  incorporated  into  a  prediction  model. 

The  four  fatality  prediction  models  developed  from  Data  Set  A  were  used 
to  pr- diet  fatality  for  Data  Set  B.  There  were  only  two  fatal  runs  in  the 


latter  set,  one  of  which  could  not  be  explained  by  any  of  the  models.  The 
other  fatality  was  accurately  predicted  by  the  two  models  which  included  AAZ 
as  a  predictor  but  not  by  the  other  two  models.  Thus,  there  is  some 
indication  that  AAZ  is  a  better  prediction  of  fatality  than  is  PSA,  but  there 
is  not  enough  evidence  for  a  firm  conclusion  to  be  drawn. 

Confidence  limits  were  given  for  the  predictions  from  each  of  the  models 
and  were  in  all  cases  fairly  wide.  In  general,  large  sample  sizes  are  needed 
to  obtain  precise  predictions  for  dichotomous  variables.  It  should  be  noted 
that  the  confidence  intervals  were  narrowest  for  the  model  based  on  AAZ  alone, 
although  the  model  based  on  AAZ  and  |PHC|  gave  slightly  better  predictions. 

It  is  not  clear  which  model  is  preferable,  but  their  predictions  are  similar 
enough  that  it  probably  makes  little  difference  which  is  chosen.  There  is 
definite  evidence  that  |  PHC|  is  of  some  value  as  a  predictor  variable,  but  its 
significance  is  small  compared  to  that  of  AAZ. 

In  order  to  obtain  the  best  possible  fatality  prediction  models,  the  two 
data  sets  were  combined  and  the  modeling  process  repeated.  In  three  of  the 
four  cases,  the  same  predictor  variables  were  selected.  The  model  based  on 
PSA  and  |PHC|,  however,  was  replaced  by  a  model  based  on  PSA  and  PHB,  where 
PHB  is  the  initial  pitch  angle  of  the  head.  Once  again,  AAZ  was  found  to  be  a 
slightly  better  predictor  than  PSA,  and  the  best  model  was  that  based  on  AAZ 
and  | PHC| . 

For  Data  Set  B,  injuries  were  classified  at  four  levels  of  severity,  with 
the  highest  level  representing  fatality.  Models  were  developed  for  each  of 
these  levels  based  on  AAZ,  which  was  found  to  be  the  single  best  predictor  of 
injury.  These  were  too  few  injuries  to  allow  the  development  of  models  with 
more  than  one  predictor  variable. 


For  each  of  the  graded  injury  models,  the  LD,.q  was  calculated.  (This  is 
the  level  of  AAZ  at  which  the  predicted  probability  is  50Z.)  Comparison  of 
the  LD^s  showed  that  the  second  and  third  severity  levels  could  not  be 
distinguished  in  this  data  set.  This  result  could  be  spurious,  since  there 
were  only  one  level-2  injury  and  two  level-3  injuries,  too  little  data  on 
which  to  base  conclusions.  However,  for  the  sake  of  this  analysis,  it  was 
necessary  to  combine  the  two  levels.  The  remaining  three  injury  prediction 
models  (slight,  moderate,  fatal)  were  plotted  and  shown  to  be  adequately 
separated  except  at  the  highest  negative  AAZ  levels. 

The  fact  that  AAZ  was  found  to  be  the  best  predictor  variable  for  each 
level  of  injury  suggests  that  the  same  mechanism(s)  may  be  responsible  for 
each  type  of  injury.  Furthermore,  the  LD,_qS  are  estimates  of  the  mean 
population  thresholds  for  each  injury  level  and  may  therefore  be  used  as  an 
initial  step  toward  the  establishment  of  a  continuous  injury  scale. 
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APPENDIX  A 


The  following  two  pages  contain  a  listing  of  the  sled  experiments  in  Data 


Set  A.  Periods  indicate  missing  values.  A  description  of  the  variables  is 
given  in  Table  1 . 
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APPENDIX  B 


The  following  two  pages  contain  a  listing  of  the  sled  experiments  in  Data 


Set  B.  Periods  indicate  missing  values.  A  description  of  the  variables  is 


given  in  Table  1. 
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